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^sj | To Victoria and Joaquin, propagators of our happiness 



C] Abstract. In this note a Galoisian approach to build propagators 

O |' through Riccati equations is presented. The main result corresponds to 

the relationship between the Galois integrability of the linear Schrodinger 
equation and the virtual solvability of the differential Galois group of its 
C^ ' associated characteristic equation. As an application we find the propa- 

gator for a generalized harmonic oscillator that has applications describ- 
ing the process of degenerate parametric amplification in quantum op- 
tics and the description the light propagation in a nonlinear anisotropic 
waveguide. 
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^i ' 1. Introduction 

£^ , The generalized harmonic oscillator has attracted considerable atten- 

tion over many years in view of their great importance to several advanced 
quantum problems including Berry's phase, quantization of mechanical of 
systems and more (see [11] and references therein). The fact that in quan- 

k> ' turn electrodynamics the electromagnetic field can be represented as a set of 

$—i . forced harmonic oscillators makes quadratic Hamiltonians of special interest 

IS [9J, IT2l I18|, 1451 150| . A method to construct explicit propagators for the 
linear Schrodinger equation with a time-dependent quadratic Hamiltonian 
based in solutions of the Riccati equation has been presented in [7], 15], and 
as particular cases the propagators for the free particle, the harmonic oscil- 
lator, Caldirola-Kanai and more can be found in a unified manner [7] \S[ [22] . 
There has also been an increase interest in the study of Picard-Vessiot 
Theory, also known as the Galois theory for ordinary linear differential equa- 
tions, where differential equations of such type are analysed throughout their 
Galoisian structure. This Galoisian structure depends on the nature of the 
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solutions of the differential equation; for instance, one obtains some kind of 
solvability (virtual) for the Galois group whenever one obtains Liouvillian 
solutions, and in this case one says that the differential equation is inte- 
grable. This means for example that when one obtains Airy functions, the 
differential equation is not integrable, while when one obtains Jacobi elliptic 
functions, the differential equation is integrable, and for instance one gets 
virtual solvability of its Galois group. 

The aim of this paper is to establish a Galoisian approach to the tech- 
niques given by Suslov et.aJL, see [Tj, [HI [SI HOj - To study Liouvillian solu- 
tions for linear second order differential equations, as well the integrability of 
their associated Riccati equations, we use Kovacic algorithm (see [21 j ) and 
an algebrization procedure (see [T], [3]) . These tools were applied to study 
differential equations incoming from physics; in particular the integrability 
analysis of the one-dimensional linear Schrodinger equation has been studied 
in pi Sl- 
As an application of this Galoisian approach we show how to find ex- 
amples of explicit propagators through Liouvillian solutions for the linear 
second order differential equation associated to a Riccati equation. As a 
first example of how to incorporate the Galoisian approach we consider the 
following (quadratic) generalized harmonic oscillator: 

id t ijj = Hip (1) 

Hit) = ( l + -cos(2ujt) )p 2 + T ^- ( 1 cos(2ut))x 2 (2) 

2m \ uj J 2 \ uj J 

+— svn{2ujt){px + xp), 

The oscillator ([I])-© might be introduced for the first time by Takahashi |47j 
in order to describe the process of degenerate parametric amplification in 
quantum optics (see also [23l [2U [26], [27J S3 SB, HZ] ) • The Hamiltonian $2§ 
had also been considered by Angelow and Trifonov [4], [5] in order to describe 
the light propagation in a nonlinear anisotropic waveguide. In [llj the 
authors constructed the propagator for ([1])- ([2]) using non-periodic solutions 
of Ince's equation (compare with classical results of periodic solutions |28j); 
they have motivated the study of non-periodic solutions of Ince's equation 
that seems to be missing in mathematical literature. In fact, they present 
a non-periodic solution that allows us to write the propagator explicitly; 
the fact of the solution being non-periodic is fundamental. However, the 
authors don't get into the details of how they found the solution. In this 
note we present in detail how to find this non-periodic solution by using a 
combination of kovacic algorithm |21j and an algebrization procedure (see pQ, 
\3\). We believe this approach can be extended to the study of propagators 
of other generalized harmonic oscillators, but here we restrict ourselves to 

©-©■ 

This note is organized in the following way: 
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• Section 2 contains a brief description of the basic theory concerning 
of construction of explicit propagators using the references O I10|, 

im © El HH HSl S2l E], anda short summar y of the Pi card- 
Vessiot theory is also presented, which was written according to 

• Section 3 contains the main results of this note; it corresponds to a 
Galoisian approach of propagators and is divided in two parts. The 
first part is devoted to a theoretical Galoisian approach to propaga- 
tors starting with Riccati and second order differential equations. 
The result given here relates the Galois integrability of the linear 
Schrodinger equation with the virtual solvability of the differential 
Galois group of its associated characteristic equation. In the second 
part is presented the Galoisian analysis of second order differential 
equations to do the construction of new propagators, as well their 
Green functions through characteristic equations and Riccati equa- 
tions. Also, we compute the differential Galois group associated to 
each propagator, which corresponds to the differential Galois group 
of the characteristic equation. 

• It has been included an appendix, written according to [Tl |2"1 [3], 
which describes the Kovacic algorithm and the Hamiltonian alge- 
brization, tools used to compute the propagators in Section 3. 

For suitability, along this paper d x denotes ^-, for higher order derivation 
(9" denotes Jpj- and by dta(0) we mean dta(t)\t=o- 

2. Theoretical Background 

2.1. Differential Galois theory. Theoretical aspects of Galoisian ap- 
proach to differential equations are summarised in this subsection which is 
based on the references [19, 20, 2TJ, \'30\, 148] , Moreover, applications of dif- 
ferential Galois theory to mathematical physics can be found in [H [21 EH I32L 
I33L 142] , The Galois theory of differential equations, also called Differential 
Galois Theory has been developed by Picard, Vessiot, Kolchin and currently 
by a lot of researchers. In particular, we focus in the Galois theory of lin- 
ear differential equations, also known as Picard- Vessiot theory. Following 
[3] and also [T], we present here an algebraic model for functions and the 
corresponding Galois theory. 

Definition 1 (Differential Fields). . Let K be a commutative field of 
characteristic zero. A derivation of K is a map d x : K — > K satisfying 
d x (a + b) = d x a + d x b and d x (ab) = d x a -b + a- d x b for all a,b G K . We then 
say that (K,d x ) (or just K, when there is no ambiguity) is a differential 
field with the derivation d x . 

We assume that K contains an element x such that d x (x) = 1. Let C 
denote the field of constants of K: 

C = {ce K\d x c = 0}. 
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It is also of characteristic zero and will be assumed to be algebraically closed. 
The coefficient field for a differential equation is defined as the smallest 
differential field containing all the coefficients of the equation. 
Due to we will mostly analyze second order linear homogeneous differential 
equations, i.e equations of the form 

C := d x y + ad x y + by = 0, a,b G K. 

so the rest of the theory will be explained in this context. 

Definition 2 (Picard-Vessiot Extension) . Let L be a differential field 
containing K (a differential extension of K). We say that L is a Picard- 
Vessiot extension of K for C if there exist two linearly independent 2/1,2/2 € L 
solutions of £ such that L = K(y\,y2) (i.e L = K(yi,y2,d x yi,d x y2)) and L 
and K has the same field of constants C. 

In what follows, we choose a Picard-Vessiot extension and the term 
"solution of £" will mean "solution of C in L". So any solution of C is 
a linear combination (over C) of y\ and y<i- 

A A-automorphism a of the Picard-Vessiot extension L is called a differ- 
ential automorphism if it leaves K fixed and commutes with the derivation. 
This means that a(d x a) = d x (a(a)) for all a € L and Va € K, o~(a) = a. 

Definition 3 (Differential Galois Groups). The group of all differential 
automorphisms of L over K is called the differential Galois group of L over 
K and is denoted by DGal(L/A'). 

Given a S DGal(L/A'), we see that {o~yi,ay2} are also solutions of C 
Hence there exists a matrix 



A * = (c d) GGL ( 2 ' C )' 



such that 

o-((yi 2/2)) = (er(yi) <j(y 2 )) = (yi 2/2) A a . 
As a commutes with the derivation, this extends naturally to an ac- 
tion on a fundamental solution matrix of the companion first order system 
associated with C. 

a ( ( J/i V2 X\ = ( o-(yi) a(y 2 ) \ = ( y\ y 2 \ A 
\\d x yi d x y 2 ) ) ~ \a{d x yx) a(d x y 2 )J ~ \d x y± d x y 2 ) a ' 

This defines a faithful representation DGal(L/A) — >• GL(2,C) and it is 
possible to consider DGal(L/A) as a subgroup of GL(2,C). It depends on 
the choice of the fundamental system {j/1,2/2}, but only up to conjugacy. 
Recall that an algebraic group G is an algebraic manifold endowed with a 
group structure. Let GL(n, C) denote, as usual, the set of invertible n x n 
matrices with entries in C (and SL(n, C) be the set of matrices with deter- 
minant equal to 1). A linear algebraic group will be a subgroup of GL(n, C) 
equipped with a structure of algebraic group. One of the fundamental results 
of the Picard-Vessiot theory is the following theorem. 
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Theorem 4. [20] The differential Galois group DG&1(L/K) is an alge- 
braic subgroup o/GL(2,C). 

In fact, the differential Galois group measures the algebraic relations 
between the solutions (and their derivatives). It is sometimes viewed as the 
object which should tell "what algebra sees of the dynamics of the solutions" . 

In an algebraic group G, the largest connected algebraic subgroup of G 
containing the identity, noted G°, is a normal subgroup of finite index. It 
is often called the connected component of the identity. If G = G° then G 
is a connected group. Furthermore, when G° satisfies some property, we say 
that G virtually satisfies this property. For example, virtually solvability of 
G means solvability of G°. 

Theorem 5 (Lie-Kolchin Theorem). Let G C GL(2,C) be a virtually 
solvable group. Then G° is triangularizable, i.e it is conjugate to a subgroup 
of upper triangular matrices. 

The following corresponds to the list of the Algebraic Subgroups o/SL(2, C), 
see [201 ED SS] for further details. 

• Reducible subgroups. These are the groups which leave a non-trivial 
subspace of V invariant. They are classified in two categories. 
— Diagonal groups: 



the identity group: 

{e} = 
the n— roots group: 



1 
1 



G [n 



c 
c" 1 



the multiplicative group: 

'c 
.0 c 



<G m = { ( n _! , c e 



Triangular groups: 
the additive group: 



(n f], 



,0 1 
the n— quasi-roots: 



v V C 

the Borel group: 

c d 



_i , ceC, deC 
c 



Irreducible subgroups 
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The infinite dihedral group (also called meta-abelian group): 

o A). - c 'HO o). "^ 

and its finite subgroups ID>2n (where c and d spans the n-th 
roots of unity). 
— Finite irreducible (primitive) groups: 
the tetrahedral group A A 2 of order 24, 
the octahedral group 5 4 2 of order 48, 
and the icosahedral group A 5 2 of order 120. 

Definition 6 (Integrability of C). We say that the linear differential 
equation C is (Liouville) integrable if the Picard- Vessiot extension L D K is 
obtained as a tower of differential fields K = Lq C L\ C • • • C L m = L such 
that Li = Lj_i(r/) for i = 1, . . . , m, where either 

(1) rj is algebraic over L.j_i, that is n satisfies a polynomial equation 
with coefficients in Lj_i. 

(2) 7/ is primitive over Li_\, that is d x rj G Lj_i. 

(3) r\ is exponential over Lj_i, that is d x n/n £ Li—\. 

We remark that the usual terminology in differential algebra for inte- 
grable equations is that the corresponding Picard- Vessiot extensions are 
called Liouvillian. The following theorem is due to Kolchin. 

Theorem 7. The equation C is integrable if and only if DG&l(L/K ) is 
virtually solvable. 

The following result is taken from [2] and corresponds to the transfor- 
mations involving Riccati equations and second order differential equations. 

PROPOSITION 8 (Riccati transformations, [2]). Let K be a differential 
field and let consider functions ao(x), a\{x), a,2(x), r{x), p(x), bo(x), b\(x) 
belonging to K that, for simplicity, will be denoted without their explicit 
dependence on x. Consider now the following forms associated to any second 
order differential equation (ode) and Riccati equation: 

(i) Second order ode (in general form): 

8 2 x y + b x d x y + b y = 0. (3) 

(ii) Second order ode (in reduced form): 

dli = PC- (4) 

(iii) Riccati equation (in general form): 

d x v = a + a\v + a 2 v 2 , a 2 / 0. (5) 

(iv) Riccati equation (in reduced form): 

d x w = r — w , (6) 
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Then, there exist transformations T, B, S and R leading some of these 
equations into the other ones, as showed in the following diagram: 

d x v = oq + a\v + a 2 v 2 »- d x w = r — w 2 



B 



R 



d 2 y + hd x y + b y = »- ti%£ = p£. 

The new independent variables are defined by means of 

d x a 2 a x \ 1 1 d x y 

2a% 2a 2 J a 2 a 2 y 

S: y = i e-y^, R: W = & -f, 

and the functions r, p, bo and b\ are given by 

r = ~5 (°o + °i a + a 2« 2 _ 9 x a) , (7) 

d x a 2 , ai \ _ 1 



V 2a 2 2a 2 J a 2 

h = - ( a\ + -2-2- ) , b = a a 2 , (9) 

P = r = ^; + — -bo. (10) 

Remark 9. From Proposition 0, it is recovered the well known result 
in differential Galois theory (see for example [48JJ: Riccati equation has an 
algebraic solution over the differential field K if and only if its associated 
second order differential equation has two independent Liouvillian solutions 
(the differential Galois group of the second order differential equation is 
virtually solvable). 

2.2. Propagators and Green Functions. In this section, as well in 
the rest of the paper, we follow [10] considering the one-dimensional time- 
dependent Schrodinger equation for an harmonic oscillator 

idtip = Hip, H = a(t)p 2 + b(t)x 2 + c(t)(px + xp), p = id x . (11) 

The Schrodinger equation (jllH can be written as 

id t i) = -a (t) d%tp + b (t) x 2 ^ - ic (t) V> - i-^-xd x ^. (12) 

We start considering the Riccati equation 

d t a + b(t)+2c(t)a + 4a(t)a 2 = 0, (13) 

where a(t), b(t) and c(t) are elements of a differential field K, with coefficient 
field C By Proposition [8] we can transform the Riccati equation ( [13]) 
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through the change of variable 

_ 1 dtfi (t) c(t) 
a[t) ~ Aa(t) //(t) 2a(tY [ } 

into the second order differential equation 

d 2 fi - r(t)dt(i + 4a(t)[i = 0. (15) 

Moreover, by Remark^ the differential Galois group of the differential equa- 
tion (|15p is virtually solvable if and only if the Riccati equation (|13p has an 
algebraic solution over the differential field K. Furthermore, by Proposition 
departing from the differential equation (|15p we can arrive at the Riccati 
equation (|13p through changes of variables. The following lemmas show how 
can we construct propagators based on explicit solutions in (fl3|) and (fI5|) . 



Lemma 10. [7 ], I10|, [44] Given a{i), bit) and c(t) piecewise continuous, 
there exists an interval I of time where the following (Riccati-type) system 

d t a + b{t) + 2c{t)a + Aa{t)a 2 = 0, (16) 

d t (3 + (c(t) + 4a(t)a(t))/3 = 0, (17) 

d a + a(t)[3 2 (t)=0, (18) 

has as a fundamental solution in terms of solutions of the following (char- 
acteristic) equation 

d 2 t ii - t (t) d t fi + 4xr(t)(i = (19) 

/ n d+a , . , c 2 c ( ' c\a d t c\ ,„„. 

rt= — "4c, a (t) = ab + - + - p M 20 

given by: 

a ° (i) -4a(t) /, (t) ~Mt)' (21) 

#>(<) = \w, (22) 

Mow 

7o(t) -2Mo)MO + Mco (23) 

provided that /j,q and Hi are standard solutions of $19\ ) -II20\ ) with /j,q (0) = 0, 
dtfio (0) = 2a (0) / 0, and m (0) / 0, S^Mi (°) = °> and S o (0) = -£o (0) = 
<7 (0) / (2a (0)), ko (0) = 0. Further, the following asymptotics hold: 

i c (o) Mo) 

ao W = M0)i ~M0J-8aW) + {l (M) 

^«-sk + Si +o «' (25) 

^ W -4^^ + 4^-8^0T + O W < 26 > 

as t — 7- 0, wi/i a(t), 6(t) and c(t) sufficiently smooth. 
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Lemma 11. [7J, I10|. I44L I45j The Green function, or Feynman's propa- 
gator, corresponding to the Schrodinger equation [77] can be obtained as 

yj = G (x, y, t) = l e i(ao(i)x 2 +/3o(t)xy+7o(t)j/ 2 ) 9 (27) 

yj2mfXo (t) 

where uq (t) , j3$ (t) and 70 (t) are solutions of the Riccati-type system. Then 
the superposition principle allows us to solve the corresponding Cauchy ini- 
tial value problem: 



ip(x,t)= G(x,y,t)ip(y,0)dy 

J — 00 

for suitable data ip(x, 0)= <p(x). Further as t — )■ 0, 



G (x, y, t) ~ — ; exp 

V ; y/2ma(0)t 



\x-y) 2 
% 4a(0)t 



x exp 



fdjam {x _ ? + ^l {x _ ? 

V8a 2 (0) lX V) + 4a(0) [X V) 



3. Galoisian approach to Propagators 

In this section we apply the Picard-Vessiot theory in the context of 
propagators. 

We are interested in Liouvillian solutions of the linear Schrodinger equa- 
tion ([lip ; that is, the so-called Liouvillian propagator (|27p is obtained 
through Liouvillian functions. In this way, we can give a Galoisian for- 
mulation for this kind of integrability. 

Definition 12. Linear Schrodinger equation ill]) is integrable in Ga- 
loisian sense (Galois integrable), when it has a Liouvillian propagator (21). 



Theorem 13 (Galoisian approach to LSE). Linear Schrodinger equation 
\11\) is Galois integrable if and only if the differential Galois group of the 
characteristic equation < f!5)) is virtually solvable. 



Proof. Consider K as the differential field of the characteristic equa- 
tion. Let us suppose (see Lemma LTUl) that fi$ is a Liouvillian solution of 
the characteristic equation; then \i\ is Liouvillian too. This implies that the 
differential Galois group of the characteristic equation is virtually solvable 
and a is an algebraic solution of the Riccati equation over K. Now, using 
Ho and \i\ in f3 and 7 such as in Lemma [101 we see that they are Liou- 
villian functions over K and therefore the propagator is Liouvillian. Thus, 
by Definition [12] we get that the Schrodinger equation is Galois integrable. 
Conversely, assuming that the Schrodinger equation is Galois integrable we 
see that it has a Liouvillian propagator that can be obtained through an 
algebraic solution of the Riccati equation over K. For instance, there exists 
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a Liouvillian solution ^o of the characteristic equation which implies that 
its differential Galois group is virtually solvable. □ 



Remark 14. In virtue of Theorem [73| we can construct many Liouvil- 
lian propagators through integrable second order differential equations over 
a differential field as well by algebraic solutions, over such differential field, 
of Riccati equations. This is the practical aim of this paper, which will be 
given in Section \3. 1\ 



3.1. New Propagators From Differential Galois Theory. In this 
section we use the Galoisian approach to LSE, Theorem [13l to general- 
ize examples introduced in [22] and find the explicit solution presented by 
Cordero-Suslov in [lOj . as well as to introduce new propagators. The start- 
ing point is the knowledge of the integrability, in the Picard-Vessiot sense, of 
the characteristic equation, or equivalently the existence of an algebraic so- 
lution over its differential field of its associated Riccati equation. To obtain 
the solutions of the characteristic equations we apply two practical tools of 
differential Galois theory: Kovacic Algorithm and Hamiltonian Algebriza- 
tion, see Appendix A. 

3.1.1. Characteristic equation <9 t 2 /i + t n fi = 0. We consider in this equa- 
tion as differential field to K = C(t). This is a generalization of the case 
n = 1, which was presented in (22] . However, the case n = 1 does not corre- 
spond to Galoisian integrability of the characteristic equation: the solutions 
are not Liouvillian due to the fact that they are Airy functions. Now, we 
study the integrability, in Galoisian sense, of this equation through Kovacic 
Algorithm. As in [42], we obtain three conditions for n to get virtual solv- 
ability of the differential Galois group. It follows the construction of the 
propagators related with these integrability conditions. 

(1) Let's consider the characteristic equation d\\x + \x = (n = 0). A 
basis of solutions is given by B = {sin t, cost}. Thus, the Picard- 
Vessiot extension is given by L = C(t,e lt ) and for instance the 
differential Galois group is DGal(L/K) = G m . By Proposition [U 
through the change of variable a = dtfi/fi, we obtain the Riccati 
equation ([13]) . where a = 1/4, 6 = 1 and c = 0. Considering 
Ho(t) = Aisint and fii(t) = A2Cost, we see that the conditions of 
Lemma [TOl are satisfied when Ai = 1/2. Furthermore, A2 must be 2 
to get W(fio, fii) = 1. In this way, by Theorem [TBI the Schrodinger 
equation is Galois integrable and its propagator is Liouvillian. 

(2) Let's consider the characteristic equation d^fi + fi/t 2 = (n = —2). 
A basis of solutions is given by B = {t m+l ,t~ m }, being m = ( — 1 db 
v5)/2. Since m £ Q, the Picard-Vessiot extension is L = C(t,t m ) 
and therefore DGal(L/K) = G m . In [JJ there is a complete study 
of the Galoisian structure of this equation. 
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(3) The characteristic equation 9ffi + fi/t 4 = (n = — 4). A basis of 
solutions is given by B = {tcos(l/t),tsin(l/t)}. Thus, the Picard- 
Vessiot extension is L = C(£, e l ' t ) and therefore DGal(L/K) = G m . 
Here, initial conditions for the characteristic equation are satisfied 
when t — > + , that is, fi (t) — > when t — > + . 

3.1.2. Characteristic equation d^fi + 2X tan(Xt)dtfJ, — 2fj, = 0. This equa- 
tion is a particular case of the Ince's equation 

o9 2Aw sin(2wt) n w 3 - 3w A 2 - (uj 2 X + A 3 ) cos (2ut) 

d;n -\ -. '-r-dtn -\ t t '-r -[i = 

uj + A cos (2 ui t) uj + A cos (2 u t) 

when A = uj. The general case corresponds to a parametric oscillator in- 
troduced in [11] . where the case A = oj = 1 was completely analyzed. To 
avoid triviality, we consider A 7^ for the characteristic equation which will 
be analysed with differential Galois theory. To do this, we apply the hamil- 
tonian algebrization process and Kovacic's algorithm on the characteristic 
equation, and so we obtain the solutions as well its differential Galois group. 
Afterward, we construct the corresponding propagator associated to this 
characteristic equation. We consider K = C(tanAi) as the differential field 
of the characteristic equation. 

We first consider the hamiltonian change of variable r = tan At obtaining 

$£(1-) + T iZ- j a r/ u(r) - (1+ 2 r2)2 Mr) = 0, (28) 

i.e., the algebraic form of the characteristic equation, being K = C(r) its 
differential field which is isomorphic to K = C(tan Xt). In order to apply the 
Kovacic algorithm we should reduce the equation (|28p through the change 
of dependent variable 

^ r) = yexp Bi TT^ s ) = TT^- 

Thus we obtain the reduced form 

^2 2r 2 + 4 2 2 

drv = ry, -=( lT 7 2 7 = iT^ + (iT^' (29) 

where the Picard-Vessiot extension is the same as in ()28|) . denoted by L, 
and for instance the differential Galois group will be the same for both 
equations due to the fact that their coefficients are rational functions in r, 
i.e., K = C(t). 

Applying Kovacic's algorithm we see that T = {i, — i,oo}, 07^ = or_j = 
oroo = 2, which implies that equation (I29f) can fall in cases 1, 2, 3 or 4 of 
the algorithm. We begin analyzing case one. 

Case 1. According to step 1, we check conditions (02) and (002) to obtain 
b{ = b-i = — \ and b^ = 2. In this way, a { = a_ i = ^-, a^ = -^ 3 -, 
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y/r i = \/r_ i = -v/Fqo = 0. By step 2 we obtain D = {1} through two dif- 
ferent options: 1 = a+ — af — oT_ i and 1 = a+ — a~ — a^; for instance 

we have two possibilities for ui: wi = ^rrj — ^r^ and 102 = 7+^ _*) that 

is, u; 2 = tt>2- Following step 3 we see that there does not exist a polynomial 
.Pi(t) = r + ao satisfying the algebraic equation (f32|) . In this way, we should 
go to the case 2 of Kovacic algorithm due to DGal(L/A') not being a sub- 
group of the Borel group. 

Case 2. According to step 1, we check conditions (02) and (002) to obtain 
hi = b-i = -\ and 600 = 2. In this way, E, = E-i = {2}, Eoo = {-4, 2, 8}. 
By step 2 we obtain D = {2} only for eoo = 8. In this way, we obtain 
6 = -^—- t + ^.. Following step 3 we obtain the polynomial Pz(t) = r 2 — 1 
satisfying the algebraic equation ([33]) . which leads us to obtain </> = ^^ + 
-^ri + z\ - In this way, solving the algebraic equation 

w 2 + (jxj + -{d T <j> + <f> 2 - 2r) = 0, 

we obtain two solutions for to: 

n r 2 -r + l „ r 2 + r + l 

CO- = 2 — - , Wj 



r 3 - r 2 + r - 1 ' + r 3 + r 2 + r + 1 ' 

In this way, due to y = e-> T , we have the general solution of (I29j) : 

y = Cie- arctanr (r - l)\A + T 2 + C 2 e arctanT (T + l)\/l + T 2 , 

for instance DG&1(L/K) = D m . Now, the general solution for equation (|28p 
is given by 

for instance the differential Galois group for the algebrized characteristic 
equation (j28l) is also the dihedral infinite group B^. Recalling that r = 
tan At, we get the general solution of the characteristic equation 

H(t) = Cie" At (sin At - cos At) + C 2 e At (sin At + cos At), 

which can also be written as 

fi(t) = (Ci + C2)(sinh At cos At + cosh At + sin At 

+(C2 — Ci) sinh At sin At + cosh At cos At, 

and its differential Galois group is also the dihedral infinite group, i.e.,DGal(L/K) 
Poo. Now, we find /J,o(t) and Hi(t) satisfying the conditions of Lemma [TOl 

l 
Ho (t) = sinh At cos A + cosh At sin At, C\ = C2 = - , 

Hi(t) = sinh At cos A + cosh At sin At, —C\ = C2 = -, 
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which for A = 1 corresponds to the solutions given in [11] . Finally, the 
propagator of the Schrodinger equation for this characteristic equation are 
given by 

G(x,y,t) = —^^^^^^=^^^^^^= (30) 

yj 2-7ri(cos At sinh At + sin At cosh At) 

c 2 ~ y 2 ) sm to smn to + 2xy — (x 2 — y 2 ) cos At cosh At 



x exp 



2i(cos At sinh At + sin At cosh At) 



Final Remarks 

This paper is an starting point to study the integrability of partial differ- 
ential equations in a more general sense through differential Galois theory. 
With this approach we studied the linear Schrodinger equation correspond- 
ing to a generalized (quadratic) harmonic oscillator obtaining similar results 
to [TTII22]. 
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Appendix A. Kovacic Algorithm and Hamiltonian Algebrization 

This appendix contains the basic information related with two practical 
tools of differential Galois theory, which have been used to compute the 
propagators in Section [3j 

A.l. Kovacic algorithm. Kovacic in 1986 (see |21j ) introduced an 
algorithm to solve the differential equation d 2 Q = r£, where r £ <C(x). 

Each case in Kovacic's algorithm is related with each one of the algebraic 
subgroups of SL(2, C) and the associated Riccatti equation 

d x v = r — v = (\/r — v) (y/r + v) , v = — r— ■ 

There are four cases in Kovacic's algorithm. Only for cases 1, 2 and 3 we 
can solve the differential equation, but for the case 4 the differential equation 
is not integrable. It is possible that Kovacic's algorithm can provide us only 
one solution (£i), so that we can obtain the second solution (£2) through 

C2 = Ci/|- (31) 
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For the differential equation given by 

d%( = r(, r = 7' s,teC[x], 

we use the following notations. 

(1) Denote by T' be the set of (finite) poles of r, F' = {c G C : t(c) = 0}. 

(2) Denote by F = F' U {oo}. 

(3) By the order of r at c G r', °(r c ), we mean the multiplicity of c as 
a pole of r. 

(4) By the order of r at oo, o (r^) , we mean the order of oo as a zero 
of r. That is o (r^) = deg(t) — deg(s). 

Case 1. In this case [y/r\ c and [y^loo means the Laurent series of -y/r at 
c and the Laurent series of y/r at oo respectively. Furthermore, we define 
e(p) as follows: if p G F, then e(j>) G {+, — }. Finally, the complex numbers 
a+, a~, a+ , a^ will be defined in the first step. If the differential equation 
has no poles it only can fall in this case. 

Step 1. Search for each c G r" and for oo the corresponding situation as 
follows: 

(co): If o (r c ) = 0, then 

(ci): If o (r c ) = 1, then 

[Vr\ e = 0, at = l. 

(02): If o (r c ) = 2, and 

r = • • • + 6(x — c)~ + • • • , then 

r n n ± l±VT+4b 

[v^J c = 0, a c = . 

(c 3 ): If o(r c ) = 2v> 4, and 
r = (a (x - c)" 1 ' + ... + d (x - c) -2 ) 2 + 6(x - c)"^ 1 ) + • • • , then 

[Vrlc = a(x - c)~ v + ... + d(x - c)~ 2 , af = -{±- + vj . 
(ooi): If o (foe) > 2, then 

[^00=°' a i= ' "^c = L 

(002): If o (foo) = 2, and r = • • • + 6x 2 + • • • , then 

(003): If o (roo) = — 2t> < 0, and 

r = (ax v + ... + df + bx 1 '" 1 + ■■■ , then 

[Vr] O0 =ax v + ... + d, and a* = -i±--vj . 
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Step 2. Find D ^ denned by 



D 



{ c£V J 



If D = 0, then we should start with the case 2. Now, if Card(.D) > 0, then 
for each n € D we search w G C(x) such that 

a; = e(oo)[^ : ] 00 + E( e ( c )[^ : ] c + «c (c) (^-cr 1 ). 
cer' 

S'tep 5. For each n £ D, search for a monic polynomial P n of degree n 
with 

<9 2 P n + 2w9 x P n + (d x u + oj 2 - r)P n = 0. (32) 

If success is achieved then £i = P n eJ ^ is a solution of the differential equa- 
tion. Else, case 1 cannot hold. 

Case 2. Search for each c£l' and for oo the corresponding situation as 
follows: 

Step 1. Search for each ccf and oo the sets E c ^$ and E^ ^ 0. For 
each c£l' and for oo we define £ c cZ and E^ C Z as follows: 

(ci): If o (r c ) = 1, then E c = {4} 

(02): If o (r c ) = 2, and r = • • • + 6(x — c)~ 2 + • • • , then 



£ c = 1 2 + AVl + 46 : Jfe = 0, ±2 J 



nz. 



(c 3 ): If o (r c ) = v > 2, then £ c = {t>} 

(ooi): If o (roo) > 2, then £«, = {0, 2, 4} 

(002): If ° (^00) = 2, and r = • • • + 6x 2 + • • • , then 

£00 = J2 + fcVl + 46 : Jfe = 0, ±2 j D Z. 

(003): If o (roo) = v < 2, then i^ = {v} 
Step 2. Find D^% defined by 

9 ( e °° ~ X/ c ) ' Ve P E ^P' P G r f • 



D=<n£ Z_i_ : ?i 

I 2 



If I? = 0, then we should start the case 3. Now, if Card(D) > 0, then for 
each n G D we search a rational function 6 defined by 






2^-^ x — c 
cer' 
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Step 3. For each n G D, search a monic polynomial P n of degree n, such 
that 

d 3 x P n + 39d 2 x P n + (3d x 6 + 3(9 2 - 4r)d x P n + {d 2 J + Wd x 6 + 6 3 - Ar6 - 2d x r) P n = 0. 

(33) 
If P n does not exist, then case 2 cannot hold. If such a polynomial is found, 
set 4> — + d x P n /P n and let w bea solution of 

w 2 + <pu + - (d x <f> + (j) 2 - 2r) = 0. 

Then £i = e^ w is a solution of the differential equation. 

Case 3. Search for each c£l' and for oo the corresponding situation as 
follows: 

Step 1. Search for each c G V and oo the sets E c / and E^ / 0. For 
each c £T' and for oo we define E C <Z7L and E 1 ^ cZas follows: 

(ci): If o (r c ) = 1, then £ c = {12} 

{pi)' If o (r c ) = 2, and r = • • • + b(x — c)~ 2 + • • • , then 

f I2fc , m 1 

£ c = ^6 + Vl + 46 : fc = 0,±l,...,±— , m = 4,6,mnZ. 

[ m 2 J 

(oo): If o (foe) = v > 2, and r = • • • + bx 2 + • • • , then 

f ^2fc m 1 

^oc = <!6+ Vl + 4& : fc = 0,±l,...,± — , TO = 4,6,12[>nZ. 

m 2 

.Step 2. Find D^0 defined by 
D = In G Z_i_ : n 



— I eoo - ^e c J ,Ve p G £ p , p G Y \ 



In this case we start with m = 4 to obtain the solution, afterwards m = 6 
and finally m = 12. If D = 0, then the differential equation is not integrable 
because it falls in the case 4. Now, if Card(-D) > 0, then for each n G D 
with its respective m, search a rational function 

~ U^x-c 
cer' 

and a polynomial S defined as 

S=H(x-c). 
cer' 
Step 3. Search for each n G D, with its respective m, a monic polynomial 
P n = P of degree n, such that its coefficients can be determined recursively 
by 

P_i = 0, P m = -P, 

P^ = -Sd x Pi - ((m - i) 9,5 - 50) Pi-(m- i) (i + 1) S 2 rP i+u 
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where i E {0,1... ,m — l,m}. If P does not exist, then the differential 
equation is not integrable because it falls in Case 4. Now, if P exists search 
co such that 

i=o v ; 

then a solution of the differential equation is given by 

where ui is solution of the previous polynomial of degree m. 

A. 2. Hamiltonian algebrization. In this section we follow PQ|3]. We 
recall that there are a lot of differential equations with coefficients that are 
not rational functions. For these differential equations it is useful, when is 
possible, to replace it by a new differential equation over the Riemann sphere 
P 1 (that is, with rational coefficients). To do this, we can use a change of 
variables. The equation over P 1 is called the algebraic form or algebrization 
of the original equation. 

Hamiltonian change of variable. A change of variable z = z(x) is 
called Hamiltonian if (z(x),d x z(x)) is a solution curve of the autonomous 
classical one degree of freedom Hamiltonian system 

f) z — f) H w 2 

az-a w H w . th H = H = w 

d x w = —o z H v ' 2 

for some V G K, where K is a differential field. Thus, z = z{x) is a 
Hamiltonian change of variable if there exists a G K such that (d x z) 2 = 
a(z). More specifically, if z = z(x) is a Hamiltonian change of variable, we 
can write d x z = y/a, which leads us to the following notation: d z = \fad z . 
We can see that d z is a derivation because satisfy d z (f + g) = d z f ' + d z g 
and the Leibnitz rules 

xtt \ at , t a a f f\ d zf ■ 9 ~ f ■ d z g 
dz(f-g) = d z f-g + f-d z g, d z - = = . 

We can notice that the chain rule is given by d z (f o g) = d g f o gd z (g) ^ 
d g f o gd z (g). The iteration of d z is given by 

d° z = i, d z = V^d z , d? = VcldJr 1 = V^d z (. . . (Vcld z )) . 

" v ' 

n times y/ad z 

We call Hamiltonian Algebrization to the algebrization process obtained by 
a Hamiltonian change of variable. 

Hamiltonian Algebrization Theorem [HE]. Consider the systems 
of linear differential equations [A] and [A] given respectively by 

8 X Y = -AY, d z Y = -AY, A = [a ij ), !=[%], Y = [y a ), Y = [y a ], 
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where a^- G K = C(z(x), d x (z(x))), a,ij G C(z) C K = C(z, y/a), 1 < i < n, 
1 < j < n, cnj{x) = aij{z{x)) and yn{x) = yn{z{x)). Suppose that L 
and L are the Picard-Vessiot extensions of [A] and [A] respectively. If the 
transformation (p is given by 

x (->• z 

^ ' 2/aOW m{z{x)) ' 
5 X i-> <9 2 

then the following statements hold: 

(1) K~K, (K,d s )~(K,d x ). 

(2) BGal(L/K) ~ DGal(L/A") C DGal(L/C(z)). 

(3) (DGal(L/i^))° ~ (DGal(L/C(z))°. 

A natural example of Hamiltonian Algebrization, and for instance of 
the introduction of the new derivative d z , is the case of second order linear 
differential equations. Consider d x y + ad x y + by = 0, using <p we obtain 
d 2 z y + ad z y + by = 0, which is equivalent to 

ad 2 z y + (jy + ^) d *y + *y = °' (34) 

where y{x) = y(z(x)), a(z(x)) = a(x) and b(z(x)) = b{x). 

In general, for y(x) = y(z{x)), the equation F(d x y, ... ,y,x) = with co- 
efficients given by ai k (x) is transformed in the equation F(d™y, . . . , y, z) = 
with coefficients given by ai k (z), where a,i k (x) = a,i k {z(x)). In particular, for 
sfa G C(z) and a,i k G C(z), the equation F(d™y, ... ,y,z) = is the Hamil- 
tonian algebrization of F(d x y, . . . , y, x) = 0. Now, if each derivation d x has 
order even, with a(z) and a,i k (z) being rational functions, then the equation 
F(d^y, . . . , y, x) = admits an Hamiltonian Algebrization for z = z(x). 
An example, that illustrate this, is given by the following linear differential 
equation: 

dl n y + a n „i(x)^ n ~ 2 y + . . . + a 2 (x)d x y + ai(x)dly + a (x)y = 0. 

Hamiltonian Algebrization Algorithm Theorem [Tj, [3]. In gen- 
eral is very difficult to find a suitable Hamiltonian change of variable, for 
this reason Hamiltonian Algebrization is a method or procedure (not an 
algorithm!). For specific families of differential equations we can obtain al- 
gorithms to apply Hamiltonian Algebrization, for example any differential 
equation 

F(d r x l y,d x b - 1 y,...,d x y,y;e^ x ,...,e x * x ) = 0, l x=p i \, Pi G Z, l<i<k, 

admits Hamiltonian Algebrization if and only if 

-p G Q*, 1 < i < k, l<j<k 
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and by means of the Hamiltonian change of variable z = e^ x we obtain the 
equation 

with coefficients in C(z). 
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